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Abstract. A detailed study of solutions to the first order partial 
differential equation H{x,y 1 z x , z y ) = 0, with special emphasis on 
the eikonal equation + z^ — h(x,y), is made near points where 
the equation becomes singular in the sense that dH — 0, in which 
case the method of characteristics does not apply. The main re- 
sults are that there is a strong lack of uniqueness of solutions near 
such points and that solutions can be less regular than both the 
function H and the initial data of the problem, but that this loss of 
regularity only occurs along a pair of curves through the singular 
point. The main tools are symplectic geometry and the Sternberg 
normal form for Hamiltonian vector fields. 



1. Introduction 
The eikonal equation in two independent variables for z = z(x,y) is 

(1.1) zl + Zy = h(x,y) 

where h is a non-negative smooth function on the plane R 2 , and sub- 
scripts denote partial derivatives. Near points (xo,yo) where h ^ all 
local solutions to this equation can be constructed by the method of 
characteristics (c/. p], Chap. 2], JTB), Chap. 10]) and questions of local 
existence, uniqueness, and regularity are fully understood. However, 
near points (xq, y®) where h vanishes the picture is much less complete. 
Our main goal is to study these questions for ( |1 . 1| ) (and more generally 
for Hamilton- Jacobi equations H(x,y, z x , z y ) = 0) in detail and show, 
contrary to some assertions in the literature, that near a point where h 
has zero as a non-degenerate local minimum value there are generally 
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infinitely many local solutions. The model case for such a problem is 

(1.2) 4 + 4 = aV + b2 y 2 

where a and b are positive. For ease of statement we discuss our general 
results in terms of this special case. We normalize the solution so 
that z(0, 0) = 0. Then (|1.2j ) has at least the four solutions z± t ± = 
l(±aa; 2 ± by 2 ). If a ^ b then (see Proposition |4.1| below) any C 3 
solution near (0, 0) is of the form 

z(x, y) = \{±ax 2 ± by 2 ) + 0((\x\ + \y\) 3 ). 

The two solutions = \{ax 2 + by 2 ) and z__ = —\{ax 2 + by 2 ) are 
unique in the sense that any C 2 solution z with leading terms z = 
\{ax 2 + by 2 ) + + \y\) 3 ) will agree with z +!+ in a neighborhood of 

(0, 0) (and on all of R 2 if z is globally defined), with a similar statement 
holding for z--. (See Q or Theorem |2.7| below for the local result and 



Oliensis [|TTJ for the global version). However, the other two solutions 



are far from unique. In Section [O] we show that for the solution 



z + _ = -l(ax 2 — by 2 ) there is an infinite dimensional family of C c 



solutions of the form z = \{ax 2 — by 2 ) + 0((\x\ + \y\) 3 ) that do not 
agree with z + - in any neighborhood of the origin (Theorem |3.2j ). (This 
gives counterexamples to some claims in the literature, cf. ill] , p. 1094], 



which would imply that z± ; ± are the only solutions.) It is also shown 
that for each k > 2 there are solutions that are C k , but not C h+1 . 
The method of proof is to use a variation of the well known method 
of characteristics for first order equations. Data is given along the two 
curves 

c + {t) = {t,t), c_(t) = (t,-t) 

(two curves rather than just one because of the nature of the singular- 
ity) and if the data is C k along these curves and agrees with the data 
for to order I < k at the origin then there is a surprising regularity 
phenomenon. The solution will be C k+l on R 2 \ {xy = 0}, as expected 
from the method of characteristics, but when considered as a function 
on all of R 2 the regularity is only C n with n := |~min{ , ^a+t 11 — 1 
(where [•] is the ceiling function) and so there is a drop of regularity 
along the coordinate axes. In general, if the data is C k then the solution 
z is C k+l off the coordinate axes, but on the coordinate axes the reg- 
ularity is determined by the order of contact of z with the "standard" 
solution z + _. 

This existence, lack of uniqueness, and jump of regularity along a pair 
of distinguished curves is not specific to the equation ( |1.2|) z 2 + z 2 = 
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a 2 x 2 + b 2 y 2 but generalizes to equations z 2 + z 2 = h(x, y) where h is of 
the form 

(1.3) h(x, y) = a 2 x 2 + b 2 y 2 + 0((\x\ + \y\) 3 ) 

(and more generally equations of the form H(x,y, z x , z y ) = 0) under 
the extra assumption that the numbers a and b are linearly independent 
over the rational numbers. The proof is based on using the Sternberg 
normal form for such an equation (and the existence of this normal 
form is only guaranteed when a and b are linearly independent over 
the rationals) to reduce the calculations to manageable proportions. 

The condition that a and b are linearly independent over the rational 
numbers is exactly the condition to insure that ( |1.3|) has a solution in 
formal power series with leading terms z = |(ax 2 — by 2 )+0((\x\ + \y\) 3 ). 
We include these calculations and use them to show that if a and b are 
linearly dependent over the rationals and m and n are integers such 
that ma — nb = and m + n > 4, then z 2 + z 2 = ax 2 + by 2 + x m y n 
has no C k solution for k > m + n with leading terms z = \{ax 2 — 
by 2 ) + 0((|x| + 1 2/ 1 ) 3 ) - Thus the independence condition on a and b is 
both necessary and sufficient for the existence of smooth saddle type 
solutions. 

A secondary goal of this paper is to advocate the use of differential 
geometric methods, especially for differential forms, symplectic geom- 
etry, and normal forms such as the Sternberg normal form, in working 
with first order equations. Thus we have included some expository 
material in Section |2| about symplectic geometry and its application to 
the method of characteristics for first order equations. 

The article can be summarized as follows. In Section |2| we present 
the basic methods of symplectic geometry as applied to the method of 
characteristics for first order equations of the form H(x,y, z x , z y ) = 
in two independent variables. For several reasons the theory is easier in 
two dimensions and our hope is that this will be useful in understanding 
the method in a concrete setting. This section also gives a proof of 
a slight generalization of a Theorem of Bruss [|j] on the existence of 
concave and convex solutions to z 2 + z 2 = h(x,y). 

In Section || we first give a detailed discussion of the model equation 
z x + z y = fl2x2 + b 2 y 2 - This leads to the existence, non- uniqueness 
and lack of regularity results already mentioned. This is followed by 
a generalization of these results to equations H(x,y, z x , z y ) = 0. The 
same type of existence, non-uniqueness and lack of regularity holds, but 
we were not able to give quite as precise an analysis of the regularity. 
The proofs make essential use of both symplectic geometry and the 
Sternberg normal form. 
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Section |] looks at solutions to z 2 + z 2 = h(x, y) from the point of 
view of formal power series. This leads to examples of equations that 
do not have any smooth solutions of saddle type. 

In an appendix we include precise statements of two of the geometric 
tools we use: the stable submanifold theorem and the Sternberg normal 
form. 

Finally we mention that equations of the type z 2 + z 2 = h(x, y) and 
more generally H(x,y, z x , z y ) = have been of interest in computer 
vision and related fields because of the "shape from shading problem" , 
where the goal is to reconstruct a surface z(x, y) given a gray-scale 
image of it. If the surface is matte then the intensity of the reflected 
light / can be modeled as being proportional to the scalar product of 
the surface normal and the illumination direction, / ~ h ■ v. If the 
viewing and the illumination directions coincide, and we choose this 
direction as the ^-direction, then n = (—z x , —z y , 1)/ ^1 + z 2 + z 2 and 
we get, after scaling, 

I(x, y)=n-z= * j , 

V 1 + z l + z v 

and equation ( [Lip can be recovered by setting h(x, y) = l/I(x, y) 2 — 1. 
Note that in applications the data function typically has quite low reg- 
ularity. In this case using viscosity solutions has proven to be effective, 
see U and 0. 

2. Application of symplectic geometry to the eikonal 

equation 

2.1. Review of the method of characteristics and construction 
of concave and convex solutions. Let R 4 have coordinates x, y, p, q. 
Then the symplectic form on R 4 is 

uj := dp A dx + dq A dy. 

Let iV 2 C R 4 be an imbedded surface. Then N 2 projects on an open 
set U C R 2 iff there are continuous functions p(x, y) and q(x, y) so that 

N 2 = {(x,y,p(x,y),q(x,y)) : (x,y) G U}. 

The submanifold iV 2 is a jet of a function iff there is an open set 
[/CR 2 and a function z G such that 

N 2 = {(x,y,z x (x,y),z y (x,y)) : (x,y) G U}. 

It is clear that if iV 2 is a jet of a function, then it projects, but the 
converse is not true. The following is a standard result, but we include 
the short proof for those not familiar with the differential geometric set 
up. 
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2.1. Proposition. Assume that N 2 C R 4 is a simply connected two 
dimensional submanifold o/R 4 of smoothness class C k for some k > 1. 
Then N 2 C R 4 zs toe jet o/ some function ( which will be unique up to an 
additive constant) if and only if it projects over some open set [/CR 2 
and the restriction of the symplectic form to N 2 vanishes. (That is, 
u(X,Y) = for all vectors tangent to N 2 .) Moreover, if N 2 is of class 
C k and is the jet of a function z, then z is of class C k+1 . 

Proof. Assume that N 2 projects over U so that N 2 = {(x, y, p(x, y), q(x, y)) 
(x,y) G U}. As N 2 is simply connected the same is true of U. Using 
coordinates on N , we see that the restriction of the symplectic 
form to iV 2 is 

uj = dp A dx + dq A dy = (p x dx + p y dy) A dx + (q x dx + q y dy) A dy 

= (Qx ~ Py) dx A dy. 

Therefore the restriction of uj to iV 2 vanishes iff p y = q x . But as U is 
simply connected this is exactly the condition that there is a function 
z (unique up to an additive constant) so that p = z x and q = z y . If N 2 
is of class C k then p(x,y) and q(x,y) are C k functions of (x,y). Thus 
z is of class C k+l . □ 

A two dimensional surface iV 2 C R 4 is a Lagrangian surface iff 
the restriction of oo to iV 2 vanishes. 

Given a function H : R 4 — > R the characteristic vector field of 
H (also called the symplectic gradient) is the unique vector yield 
on R 4 such that for all vectors X 

uj(£ H ,X) = -dH(X). 

Then a simple calculation yields that 

, n , s , dH d OH d dH d OH d 

(2-1) Zh = ^^ + 



dp dx dq dy dx dp dy dq 
An immediate consequence of the definition of is that 

dH(Z H ) = -u(Z H ,Z H ) = 

and so H is constant on the integral curves of 

The following is the differential geometric justification of the method 
of characteristics for the Hamilton- Jacobi equation H(x,y, z x , z y ) = 0. 
Given a smooth function H: R 4 — > R and P G R 4 , the restriction of 
dH to T(R 4 )p is denoted dHp. 

2.2. Proposition. Let N 2 be a simply connected two dimensional sub- 
manifold of R 4 . We also assume that the following non- degeneracy 
condition holds: 
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(ND) The set of points P G iV 2 where dHp ^ is dense in N 2 (here 
dHp is being viewed as a linear functional on R 4 and not just 
restricted to T(N 2 ) ). 

Then N 2 is the jet of a solution to the equation H(x,y, z x , z y ) = if 
and only if the following three conditions hold 

1. N 2 C{P: H{P) = 0}, 

2. iV 2 projects over some open set U C R 2 ; and 

3. the characteristic vector field is tangent to N 2 at every point 
ofN 2 . 

If N 2 satisfies these conditions and is a C k submanifold, then the so- 
lution z is of class C k+1 . 



2.3. Remark. For our future applications it is important to realize that 
although if is a smooth function, the zero set {H = 0} need not be 
a smooth submanifold of R 4 . For a, b > we will be interested in 
H(x, y,p, q) = p 2 + q 2 — a 2 x 2 — b 2 y 2 . The zero set {H = 0} is then a 
cone that is singular at (0, 0, 0, 0). However, if N 2 is a smooth surface in 
{H = 0} that is everywhere tangent to and which projects onto an 
open set in R 2 , then iV 2 is the jet of a solution to z 2 + z 2 = a 2 x 2 + b 2 y 2 . 
This is because dH only vanishes at one point and thus the set of points 
on TV 2 where dH ^ is dense. □ 

It is useful to give a name to submanifolds that satisfy two of the 
conditions of the proposition: 

2.4. Definition. Let H: R 4 — > R be a C k function with k > 2 and 
characteristic vector field Then a connected two dimensional sub- 
manifold iV 2 of class C 1 is an invariant Lagrangian surface iff 

1. iV 2 is a Lagrangian surface (i.e. the restriction of u to iV 2 van- 
ishes), and 

2. The characteristic vector field is tangent to iV 2 at all points of 
N 2 . □ 



Proof of Proposition \2.£\ . First assume that the three conditions hold. 
Then in light of Proposition |2TT| it is enough to show that the restriction 



of u to iV 2 is zero. If we are at a point P G iV 2 where dH ^ 0, then 
from the definition of or the formula ( |2.1|) it follows that 
Let X be a tangent vector to iV 2 at P linearly independent from £h{P)- 
Then {£h(P),X} is a basis for the tangent space T(N 2 ) P . Using the 
definition of and that dH(X) = (as H\ N2 = 0) 



u(£ H ,X) = -dH(X)=0. 
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Thus u restricted to N 2 vanishes at P. But we are assuming that the 
set of points P where dH does not vanish is dense, thus by continuity 
the restriction of to to iV 2 is zero on all of N 2 . 

Conversely, if N 2 is the jet of a solution, then clearly iV 2 C {P : 
H(P) = 0} and N 2 projects onto an open subset of R 2 and also H\ a = 
0. Let P £ N 2 . Then lu(X,Y) = for all X,Y e T(N 2 ) P . But 
a calculation shows that if Z is any vector tangent to R 4 at P with 
u(Z,X) = for all X £ T(N 2 ) P , then Z £ T(N 2 ) P . But for X £ 
T(N 2 )p we have dH(X) = so (as above) u(£ H ,X) = -dH(X) = 0. 
Thus £h £ T(N 2 ) P . This completes the proof. □ 

This result makes the geometry of the method of characteristics 
clearer than the classical presentations. Let H be of class C k for k > 2. 
Then the formula (|2.1|) makes it clear that the vector field £h is of 
class C* -1 . Let $f be the flow of £ H . That is, $f (P) = P and 
c{t) := &^(P) is an integral curve of the vector field Then (see p|, 
Thm 1, p. 80] or 0, p. 230]) the map (t, P) ^ $f (P) is C7 fc_1 . Now let 
c: (a, 6) — > R 4 be a curve of class C z (/ > 1) and let n: R 4 — > R 2 be 
the projection 7i(x,y,p,q) = (x,y). Assume 

1. H(c(s)) = 0, 

2. c projects over an imbedded[] curve of R 2 (that is s i — > tt(c(s)) is 
an imbedded curve in R 2 ), and 

3. at all points c(s) the vectors 7r*c'(s) and ti*£,h{c(s)) are linearly 
independent. 

Then let P(s,t) := $^(c(s)). By the implicit function theorem for r 
small enough the submanifold iV 2 := {F(s,t) : (s,t) £ (a, 6) x (— r, r)} 
will project over an open subset U of R 2 . Also, as H is constant 
along the flow of and H(c(s)) = 0, we have that H(F(s,t)) = 
H(&j*(c(s))) = 0. Thus PJjya = 0- By construction ^ is tangent 
to iV 2 and so iV 2 is the jet of a solution z to H(x, y, z x , z y ) = by 
Proposition |2.2|. The regularity of A^ 2 is (J mm i k - 1 >i} an( j therefore the 



regularity of 2 is (7 min { fc . / + 1 }. 

2.2. Solutions near critical points of H. We now look at the more 
interesting case of finding solutions near a critical point, P , of H . As- 
sume that dHp = and by adding a constant to H we can assume that 
H(Pq) = 0. Then near Po the set {H = 0} need not be a submanifold 
of R 4 . Assume that Pq is a non-degenerate critical point so that the 



: A curve c is imbedded iff it is the image of a smooth injective map 7 : / — > R 2 , 
with / an interval in R, so that the velocity vector 7' never vanishes and so that 
the topology of as a subset of R 2 is the same as the topology induced by 7 ( 
i.e. 7 is a homeomorphism) . 
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Hessian of H at P$ is non-singular. To make the notation easier we 
assume that P — (0,0,0,0). Then using the form Q2.ip we see that 
the linearization of the characteristic system for this vector field at the 



origin 


is 




(2.2) 








~x(t) 




d 


yit) 




dt 


P (t) 











xp 



H, 
H, 



H 



xq 



y<-i 



H, 



pp 



pi 



H, 



pq 
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-H 
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xy 



H 



yy 



-H 

-H, 



xp 



IIP 



-H, 
-H, 



x:q 



y<-i 





x(t) 




x(t) 




yit) 


= L 


yit) 




pit) 




pit) 




_q(t)_ 




q(t) 



where all the second partial derivative are evaluated at (0, 0, 0, 0) and 
this equation defines L. Let det (_D 2 if) be the determinant of the Hes- 
sian at (0, 0, 0, 0) and let 

r2 



C'2 



2H x y Hpq 2,H X q Hyp 



+ HyyHqq - H yq + H x 



The characteristic polynomial of L is then, using that det(L) 
(which can be seen by noting that L = J D 2 H, where J = 

that det(L) = det(J) det(D 2 H) = det(D 2 H)), 

det(AJ - L) = A 4 + c 2 A 2 + det(D 2 H). 
Therefore the eigenvalues are 
(2.3) 



TT TT2 

n VV n xp- 

= det(D 2 H), 
o i 2 
-h o 



so 



-c 2 



'c\ -Adet(D 2 H) 



±1 



-C 2 - V C 2 



Adet(D 2 H) 



Assuming that there is no eigenvalue with zero real part we see that 
there are exactly two eigenvaules with positive real part and two with 
negative real part. Then let be the local stable manifold for the crit- 
ical point at (0, 0, 0, 0) and N 2 the local unstable submanifold. (Loosely 
Nl is the set of points P G R 4 so that lim^ ±00 $f(P) = (0,0,0,0). 
See Appendix A.l| .) Because of the condition on the eigenvalues of 
L both and are two dimensional. If H is of class C k for 
k > 2 then the flow $f is of class C k ~ l . Therefore iV 2 and N 2 
are C k ~ 1 submanifolds of R 4 (see [^, Thm 5.20, p. 49]). Then for 
a point P G N± we have, by the invariance of H under the flow, that 
H{P) = lim t _ ±00 H($f{P)) = H(0, 0, 0, 0) = 0. Thus H\ n2 = 0. Also 

it is clear that the characteristic vector field £h is tangent to N±. The 
condition of Proposition ^]2| that does not hold automatically for N± 
is that of being locally projectable. Summarizing: 

2.5. Theorem. Let H : R 4 — >• R be a C k function with k > 2. As- 
sume that H has a non- degenerate critical point at (0, 0, 0, 0) with 
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if (0, 0, 0, 0) = and assume that no eigenvalue of L as defined above 
has zero real part. Then there are exactly two eigenvalues Ai, A2 of L 
that have positive real part. Let e\ and e 2 be the eigenvectors corre- 
sponding to Ai and A 2; let Af be the subspace spanned by e\ and e 2 and 
let 7r: R 4 — > R 2 be the projection n(x,y,p,q) = (x,y). Assume that 
the restriction vr]^: J\f — > R 2 is nonsingular. Then near (0,0,0,0) 
the unstable submanifold N 2 of is the jet of a C solution z to 



II {x. y, z x . Zy 
manifold 



1 = 0. An analogous statement holds for the stable sub- 

oftH- 



Proof. The tangent space to N 2 at (0, 0, 0, 0) is M so if tc\ m : M —> R 2 
is nonsingular the implicit function theorem implies that near (0, 0, 0, 0) 
the submanifold N 2 will project onto an open subset of R 2 . As noted 
all the other hypothesis of Proposition 22 hold. This completes the 
proof. □ 

As an example of this assume that H(x,y,p,q) = f(p,q) — h(x,y) 
where / and h have critical points at (0, 0) and the Hessians of / and 
h at (0, 0) are both positive definite. (This is a case that comes up in 
the shape from shading problem.) Then the linear map L is 













fpp fpq 
fpq fqq 



•xy 



h 



"xy 
yy 













A 
B 



where A and B are positive definite 2x2 matrices. As A is positive 
definite it has a square root Az . The matrix A2BA2 will also be 
positive definite and will consequently have positive eigenvalues. Let 
a 2 and b 2 be the eigenvalues of A2BA2 , where a, b > 0, and e 1; e 2 the 
corresponding eigenvalues, 

A^BA^e 1 



a 2 e ly 



A2BA2e 2 = b 2 e 2) 



and define 



Vl 



aA~2e x 



v 2 :-- 



A2& 2 

bA 



2e 2 



v 3 



A2e 1 
-aA 



Vl 



A2e 2 
-bA 



1 

2e 2 



By direct calculation, we find 



Lvi = ovi, Lv 2 = 6v 2 , LV3 = 



— av3, LV4 = — 6V4. 
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Thus the eigenvalues of L are ±et and ±6. The eigenvectors correspond- 
ing to the two eigenvalues with positive real part are Vi and v 2 . For 
future reference we record some elmentary facts about these eigenvalues 
and eigenvectors. The proof is left to the reader. 

2.6. Lemma. If a ^ b then there are exactly six two dimensional sub- 
spaces of the tangent space T(R 4 ) that are invariant under L, corre- 
sponding to the six ways of choosing two of the eigenvectors v 1; . . . , v 4 . 
As u(vi,Vs) 7^ and co>(v 2 ,v 4 ) ^ the subspaces span(vi,v 3 ) and 
span(v2, V4) can never be tangent to jets of functions and therefore can 
be disregarded from our considerations. Each of the remaining four 
pairs {vi,v 2 }, {v 3 ,V4}, {vi,v 4 }, and {v 2 ,v 3 } spans a two dimen- 
sional subspace of the tangent space T(R 4 ) on which uj vanishes. If 
{v.j, Vj} is any one of these four pairs and ir: R 4 — > R 2 is the pro- 
jection n(x,y,p,q) = (x,y) then 7r*Vj and ir^Vj are linearly indepen- 
dent. Therefore, if N 2 is any two dimensional submanifold tangent to 
span(vj,Vj) then, by the implicit function theorem, N 2 locally projects 
over some open neighborhood U of (0,0) in R 2 . Finally, span(vi, v 2 ) 
is the tangent space at the origin to the unstable submanifold of and 
span(v3, V4) is the tangent space at the origin to the stable submanifold 

0f£ H . □ 

The following gives a generalization of a theorem of Bruss M, who 
considered the special case of u 2 + u 2 = h(x,y). The proof uses the 
same geometric idea as Bruss's proof. 

2.7. Theorem (A. Bruss [|J). Let f(p,q) and h(x,y) be functions of 
class C k , with k > 2, and assume that f and h both have critical points 
with positive definite Hessians at (0,0). Then in the class of functions 
that are concave near (0,0) there is a unique (up to an additive con- 
stant) solution to f(z x , z y ) = h(x, y). This solution is C k and has as its 
jet near (0, 0) the stable manifold of the characteristic vector field £h 
of H(x,y,p,q) = f(p,q) — h(x,y). (Likewise, there is a unique convex 
solution, it is C k and has as jet near (0, 0) the unstable submanifold of 
in-) 

Proof. Follows from the discussion above. □ 

3. Construction of saddle type solutions near a regular 

critical point 

3.1. Construction of all solutions to z 2 + z 2 = a 2 x 2 + b 2 y 2 . We 

will now construct all "saddle point" solutions of the equation 

2 1 2 2 2 1 7,2 2 
z x + z y = a x +b y , 
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where a, b > 0. The analysis here is a model for the more general 
setup covered in Section |3]2|. This special case is also of interest as 
it is possible to be somewhat more precise about the regularity of the 
solutions. What makes this equation especially easy to analyze is that 
the components of the characteristic vector field are linear so that find- 
ing the flow involves no more than linear algebra. We will assume 
that we have a solution defined near (0, 0) with z(0, 0) = 0. If a ^ b, 
then by a formal power series argument as in the proof of Proposi- 
tion |4.i| below, we have that if z G C 3 then near (0, 0) the first few 
terms of its Taylor series are z = |(±ax 2 ± by 2 ) + 0((\x\ + \y\) 3 ) and 
if a = b we can bring z to this form by a rotation of the axes. If 



z = \(ax 2 + by 2 ) + 0((\x\ + \y 
the function is convex. Hence z 
Likewise, if z = — % 



then Theorem |2.7| applies, since 
ax 2 + by 2 ) is the unique solution. 



ax 2 + by 2 )+0((\x\ + \y\) 3 
Next we look for a solution of the form 
1 
2 

We will use the Hamiltonian H = |(jo 2 
acteristic vector field is then = 



then z 



ax' 



+ by 2 



z = -(ax 2 -by 2 ) + 0((\x\ + \y\) 3 ). 



Q 



a 2 x 2 



b 2 y 2 ) 



The char- 

* . , ' ; , \X x ~7\ I b xi . J. he 

ox oy op a oq 

integral curves of this vector field satisfy the differential equation, cf. 



X 




"o 





1 


cf 




X 




X 


y 













1 




y 


= L 


y 


p 




a 2 













V 


p 


3. 







b 2 

















d_ 

dt 



The eigenvectors of the matrix L are 



(3.1) 



which satisfy 





1 









1 














1 









1 


Vl = 


a 


, v 2 = 







—a 


, v 4 = 












-b 









b 



avi, Lv 2 



-av 3 , 



Lva = b\A 



The jet of the function z = |(ax 2 — by 2 ) is the span of vi and v 2 . 

To construct an invariant Lagrangian surface through (0, 0, 0, 0) we 
start with two curves j± : R — > R 4 and generate the surface by moving 
these curves by the flow of the characteristic vector field. Because of the 
nature of the singularity, two curves rather than just one are required. 
Toward this end let (p + , y?„ : R — ► R be two functions so that 

1. <p± are C k functions for some k > 1 and 
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2. tp± vanish to order I at the origin for some I < k (where / need not 
be an integer). Specifically, this means that there are functions 
Tp±: R — > R so that 

(3.2) ip±(s)=lp ± {s)s l 

where Tp ± {s) are C k functions on R \ {0} and the derivatives 
j^Tp±(s) are bounded on [-1, 1] \ {0} for < j < I. 
Define two curves 7+, 7_ : R — > R 4 by 

a 2 

7±(s) = svi ± sv 2 + (p±(s)v 3 =F ^-</?±(s)v 4 . 

Then 

ff( 7± (a))=0 

for all s. Let e* L be the exponential of L so that for P e R 4 the integral 
curve of the characteristic vector field through P is £ 1— > e' L P. Define 
F± : R 2 -> R 4 by 

(3.3) 

F±(s,t) := e' L 7 ±(s) = se at Vl ± Se - 6 V 2 + ^ ± (s)e- at v 3 =F ^±(s)e w v 4 . 

We now consider F+(s, t) for s > 0. Do the change of variables (s, £) 1— > 
(w, f ) with -u, f > given by 

■u = se 1 t = In 



v = se bt 



a + b v 

S = u b/(a+b) v a/(a+b) 



In these coordinates we have that the image of F + (s, t) with s > can 
be parameterized by 



G(u, v) —uvi + vv 2 



2 

b a , a ti -b , b 



+ , U a + i 't> a + b (/9 + (-U a + b f a + b )v3 — — U a + b V a + b lf + (u a + b V a + f >)v 4 . 

3.1. Lemma. Lei </? 5e a function on [0, 00) of class C k that vanishes 
to order I < k at t — and Zet 1 > a, (3 > 0. Then the function 



E(u,v) :-- 



u a l v l3 (p{u a v f5 ), u,v > 
0, otherwise 



is C n on R 2 for all n < min{(/ + l)a — 1,(1 + l)/3}. Conversely, 
if E(u,v) is of class C n , then ip vanishes of order I at t = for all 
I > max{(n + l)/ct — 1, n/(3 — 1}. 
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Proof. The hypothesis on ip implies that ip(t) = t l ip(t) where ip has 
continuous derivatives up to order k on the open interval (0, oo) and 
is bounded on the closed interval [0, oo). Thus E can be rewritten as 
E(u,v) = u^ l+1 ^ a ~ 1 v ( ' l+1 ^ l3 ijj(u a v 13 ) for u, v > 0. The result now follows 
by direct calculation of the derivatives. □ 

This directly implies: 

3.2. Theorem. Let (p±: [0, oo) —>■ R be C k functions that vanish to 
order I < k atO. Let F± be defined by equation (\3.3[) and let N 2 (ip +1 ipJ) 
be the closure of the union of the images of F + and F_, that is, 

N 2 (tp + , <y2_) := image(F + ) U image(F_). 

Set n : = [min{ ^^ a , ^ a f f 1 ^ 6 }] — 1, where [■] is the ceiling function. 
Then N 2 (ip + ,Lp_) is a C n ~~ l submanifold o/R 4 that is the jet of a C n 
solution z(x,y) to the equation z 2 + z 2 = a 2 x 2 + b 2 y 2 . Conversely, 
any C n solution of this equation is of this form for unique ip+ and ip_ . 
Regardless of the value of I, the function z will be of class C k+1 on 
R 2 \ {xy = 0}. Thus there is a decrease in regularity of solutions from 
C k+1 to C n along the coordinate axes. 

Proof. All but the statements about the decrease in regularity along 
the coordinate axes follow from Lemma |3.1| . The two curves where 
the regularity of the surface N 2 ((p + , yj_) drops are along the lines u h- > 
G(u, 0) (which projects down onto the x-axis) and v i— > G(0, v) (which 
projects down onto the y-axis). All other points of N 2 ((p + , cp_) are C k . 
As in Proposition p.2|, this implies that z is C k+l off of the coordinate 



axes. □ 

3.2. Solutions to H(x, y, z x , z y ) = near a critical point of H. We 

now consider an equation H (x, y, z x , z y ) = near a critical point of H. 
We assume that H G H(0, 0, 0, 0) = and dH = at the origin. 
Let L be the linearization at the origin of the characteristic system 
defined by (|2.2j ) and assume that the eigenvalues of L are a, b, —a, —b 
where a, b > 0. (By ( [2.3|) the eigenvalues are of this form if they are 
real and nonzero.) Let ei, e^, e^, be the eigenvectors for a, b, —a, —b 
respectively. We make the two assumptions 

(3.4) a and b are linearly independent over the rational numbers, 
and 

{7r*ei,7r*e 2 }, {7T*e 3 , 7r*e 4 }, {7r*ei, 7r*e 4 }, {-K*e 2 , 7r*e 3 } 
are each linearly independent sets. 

(Where, as usual, 7r: R 4 — > R 2 is ir(x,y,p,q) = (x,y).) This and the 
implicit function theorem imply that, if iV 2 is any C l surface in R 4 that 



(3.5) 
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is tangent to both e\ and e 2 (or one of the other pairs {e 3 , e 4 }, {e 4 , e 4 } 
or {e 2 , 63}), then locally iV 2 projects over an open neighborhood of the 
origin in R 2 . 

The number theoretic assumption Q3.4|) is what is needed to in- 
voke the Sternberg Normal Form (Theorem |A.2| in the Appendix) 



and conclude there are local coordinates x,y,p,q centered at the ori- 
gin of R 4 with uj = dp A dx + dq A dy and a function f(u,v) with 
/(0,0) = 0, / u (0,0)_=_ /„(0,0) = 1 and H = \f{p 2 - a 2 x 2 ,q 2 - b 2 y 2 ). 
In the coordinates x, y, p, q the eigen-directions of L are given by the 
vectors v 1 ,V2,v 3 ,v 4 of ( |3.9| ) below. So in these coordinates we can 
take e\ = v 1; e 2 = v 4 , e 3 = v 3 and e 4 = v 2 . We then can com- 
pute that cj(ei,e 2 ) = cu(e 3 ,e 4 ) = u;(ei,e 4 ) = u;(e 2 ,e 3 ) = 0, but that 
w(ei, e 3 ), co>(e 2 , e 4 ) 7^ 0. This leads to a general existence result. 

3.3. Theorem. Let the origin 0/R 4 be a non- degenerate critical point 
of H such that the linearization L of the characteristic vector field £h 
at the origin satisfies the conditions t \3.J\) and (\3.b\ ). Then, with the 
notation above, each of the pairs {ei,e 2 } ; {e 3 ,e 4 }, {e 4 ,e 4 } and{e 2 ,e 3 } 
is tangent to the jet of a C°° solution to H(x, y, z x , z y ) = 0. 

Proof. Consider the four submanifolds N± ± of R 4 defined locally near 
the origin in the coordinates x, y, p, q by 

N± t± := {(x,y, ±ax, ±by) : (x,y) G U}, 

where U is a small neighborhood of (0, 0) in R 2 . On N± t ± the relations 
p = ±ax and q = ±by hold. Thus on N±± we have H = ^f(p 2 — 
a 2 x 2 ,q 2 - b 2 y 2 ) = |/(0, 0) = and therefore N^ ± C {H = 0}. A 
direct calculation using the form of the characteristic vector field in 
the coordinates x, y, p, q (see ( |3.7| ) below) shows that is tangent to 
Nl ± . Thus Nl ± is an invariant Lagrangian surface. We also have 

T(Af 2 i+ ) = span(ei, e 2 ), T(N 2 _) = span(e 3 , e 4 ), 
T(A^ 2 _) = span(ei, e 4 ), T(A^ 2 ,) = span(e 2 , e 3 ). 



Therefore the assumption ( |3.5|) implies that locally near (0, 0, 0, 0) each 
of the surfaces N± ± projects onto an open neighborhood of (0, 0) in 
R 2 . Consequently Proposition ^J] implies that N± ± is the jet of a C°° 
solution to H(x, y, z x , z y ) =0. □ 

If A^ 2 C R 4 is the jet of a solution to H(x, y, z x , z y ) = that passes 
through (0,0,0,0), then A^ 2 is an invariant Lagrangian surface in R 4 . 
The invariance of A^ 2 under the flow of the characteristic vector field 
implies that T(N 2 )o is invariant under the linearization L of This, 
coupled with the fact that to vanishes on T(N 2 ) , implies that T(N 2 ) 
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is one of the four subspaces span(ei,e2), span(e3,ezt), span(ei,e4), or 
span(e 2 ,e 3 ) of T(R 4 ) . 

The subspace span(ei, e 2 ) is tangent to + at the origin. This is the 
unstable submanifold of at (0, 0, 0, 0) and is unique in the sense that 
if N 2 is a C 1 invariant submanifold for with T(N 2 ) = span(ei, e 2 ) 
then N 2 = N 2 + in a neighborhood of the origin. To see this note that 
the restriction Ch\^ 2 nas a source at the origin (as the eigenvalues of the 
restriction L\~ 2 are a, b > 0) and so for all points P of N 2 sufficiently 
near the origin lirrit_>_ QO 3>t(-P) = (0,0,0,0). But this is exactly the 
condition that P be on the unstable submanifold. Thus N 2 = , as 
claimed. There is a corresponding uniqueness statement for the stable 
submanifold iV 2 _ . 

Call a solution z to H(x, y, z x , z y ) = with ^(0, 0) = z y (0, 0) = a 
solution with unstable jet iff its jet near the origin is the unstable 
submanifold + for Likewise, z is a solution with stable jet 
iff its jet near the origin is the stable submanifold N 2 for Our 
discussion leads at once to a uniqueness result. 



3.4. Theorem. If the hypotheses of Theorem \3.3j hold and if z is a 
solution to H(x,y, z x , z y ) = with z(0, 0) = 2^(0,0) = z y (0,0) = 
of either stable or unstable type, then z is unique in the sense that 
if z is any C 2 solution with 5(0,0) = z x (0, 0) = z y (0,0) = 0, such 
that the jet of z is tangent to the jet of z at (0, 0, 0, 0), then z = z in 
a neighborhood of the origin. (The tangency condition of the jets is 
equivalent to z xx (0,0) = z xx (0,0), z xy (0,0) = 2^(0,0) and z yy (0,0) = 
z yy (0,0).) □ 

Call a solution to H(x,y, z x , z y ) = with -2(0,0) = ^(0,0) = 
z y (0, 0) = whose jet at the origin is not the stable or unstable 
submanifold of the characteristic vector field a saddle solution. 
The motivation of this terminology is that for the eikonal equation 
z 2 + z 2 = a 2 x 2 + b 2 y 2 with a, b > the solution with unstable jet is the 
convex solution z = |(ax 2 + by 2 ) and the solution with stable jet is the 
concave solution z = —\{ax 2 + by 2 ). The other two obvious solutions 



z = ±|(ax 2 — by 2 ) (which are not unique by Theorem |3.2|) have saddle 
points at the origin. 

We now wish to investigate in detail the lack of uniqueness of saddle 
solutions. To do this it is more convenient to work in the coordinates 
x, y, p, q. To simplify notation we drop the bars and just write x, y, p, q 
but keep in mind that these are not the original symplectic coordinates 
on R 4 and that once we have constructed the jets of solutions we have 
to translate these results back to statements in the original symplectic 
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coordinates. With this in mind, we assume that H has form 

1 



(3.6) 
where 

and 



H 



2„.2\ 



£f 2 2 2 2 1.2 

f{p -a x ,q -by 



/(0,0)=0, / u (0,0) = /„(0,0) 



io = dp A dx + dq A dy. 
Then the characteristic vector field is 



JuPtt + fvQa- + a JuXtt- + o f v y-K-- 
ox ay op oq 



(3-7) i H 

The integral curves of this vector field satisfy the differential equation 
(3.8) 





X 










fu 







X 


d 


y 













fv 




y 


dt 


p 




fuO 2 













p 




q 







f v b 2 










Q 



This is not a linear system as the functions /„ and f v depend on x, y, 
p and q. However we will show that because of its special structure it 
can be treated almost as if it were linear. Letting L be the matrix of 
this system we find that its eigenvectors are, cf. (|3.1|), 



(3.9) 





1 









1 














1 









1 


Vi = 


a 


, v 2 = 







—a 


, v 4 = 












-b 









b 



and that 

Lvi = / u avi, 



-f v bv 2 , Lv 3 



-f u av 3 , Lv 4 = f v bw A . 



Thus L has a basis of eigenvectors that are independent of the variables 
x, y, p and q. For use in the statement and proof of Theorem p.9| we 
define the two curves 



(3.10) 



a 1 (t)=tv 1 , a. 2 (t) = tv 



2- 



which are just the first two coordinate axes for the basis v 1; v 2 , v 3 , v 4 . 
What is special about these curves is that while we will construct many 
surfaces that are jets of solutions and tangent to span(v 1 ,v 2 ), these 
curves will lie on all of these surfaces. 
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Let wi, w 2 , W3, W4 be coordinates on R 4 such that 



IUiVi + U> 2 V 2 + W3V3 + W4V4. 



Then 



(3.11) 



1 1 

w 1 = -(x + p/a), w 2 = - (y - q/b) , 

w 3 = ^(x-p/a), w 4 = -(y + q/b) . 

In these coordinates the differential equations for the characteristics 
become 

w[(t) = f u awi(t), w' 2 (t) = -f v bw 2 (t), 

Ws(*) = -f u aw 3 (t), w' 4 (t) = f v bw' 4 (t). 
Recall that f u (0, 0) = f v (0,0) = 1. Therefore given e > there is 
a 5 > so that if B(5) := {(u>i, w 2 , w 3 , W4) : Yli=i w i — * s the 
closed ball of radius 5 at the origin, then \f u — 1|, \f v — 1| <e < 1/2 
in B(5). We now use that if a function satisfies a differential 
equation = ch(t,w(t))w(t) on an interval [—«,«;], where c is a 

constant and \h(t, w(t)) — l\ < e, then on [—k, k], w(t) = w(0)e c6, W* for a 
function that satisfies |1 — 6\ < e. The differential equations for the 
characteristics are all of this form, so if (wi(0), w 2 (0), w 3 (0), w 4 (0)) G 
B(5) and t is so that (u>i(t), u> 2 (t), w 3 (r), w 4 (r)) G 5(5) for r between 
and t, then 

= Wl (0)e adlt , w 2 (t) = w 2 (0)e- be ^, 
w 3 (t) = w 3 {0)e- ae ^, w 4 (t) = w 4 {0)e b6it , 



where 



\l-6i\<e<- for i = 1,2,3,4. 



2 

Let ip + , if - : R — > R be two functions so that 

1. ip± are C fc functions for some k > 1 and 

2. </? ± vanishe to order / at the origin for some I < k. 

3.5. Lemma. There are unique functions ip± defined in a neighbor- 
hood of so that if ~f± are the curves 

7±(s) := svi ± sv 2 + v?±(s)v 3 + ^±(s)v 4 , 
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then 



H( 1± (s)) = 0. 



Moreover, 

1. ip± are C k functions and 

2. ip± vanish to the same order I < k at the origin that tp± do. As 
above, this means there are ip ± : R — ► R so that 



where ip±(s) are C k functions on R \ {0} and the derivatives 
■^j4'±(s) are bounded on [— 1, 1] \ {0} for < j < I. 

(In the caseH = \{p 2 +q 2 —a 2 x 2 —b 2 y 2 ) we haveip±(s) = =F(a 2 /b 2 )(p±(s) .) 

Proof. We will show how to find the derivation for tp_ being simi- 
lar. In terms of the coordinates W\, w 2 , w 3 , w 4 we are looking for ip+(s) 
so that if Wi = s, w 2 = s, w 3 = <p+(s), and w 4 = ?p+(s), then 
7+(s) = it?iVi + w 2 v 2 + W3V3 + W4V4 is the required curve. Writ- 
ing H in terms of the coordinates x,y,p,q (related to the coordinates 
Wi, w 2 , W3, W4 by (|3.11|) ) we have H = jf( y p 2 — a 2 x 2 ,q 2 — b 2 y 2 ). Then we 
want ip+ so that H(^ + (s)) = ^f(—Aa 2 sip + (s), —4b 2 stfj + (s)) = 0. Recall- 
ing the assumptions on / (i.e. /(0,0) = and / u (0, 0) = f v (0,0) = 1) 
we have by the implicit function theorem that there is a smooth func- 
tion v(-) defined in a neighborhood U of with v(0) = 0, v'(0) = — 1 
so that f(u,v(u)) = 0. Then near (0,0) we have {(u,v) : f(u,v) = 
0} = {(u,v(u)) : u G U}. Therefore near s = we are trying to 
solve — 4b 2 sip + (s) = v(— 4a 2 s<^ + (s)) for ip + (s). As v(u) vanishes at 
the origin it can be expressed as v(u) = uv(u), where v is a smooth 
function and v'(0) = —1 implies v(0) = — 1. Therefore i[>+(s) = 
(a 2 /b 2 )(p + (s)v(— 4a 2 sip + (s)). Note that if f+(s) vanishes to order I 
at s = 0, then this formula makes it clear that the same is true for 
ip+(s). Finally, if H = |(p 2 + q 2 — a 2 x 2 — b 2 y 2 ) then f(u,v) = u + v and 
so v(u) = —u, v(u) = —1 which implies i/)+(s) = —(a 2 /b 2 ){p + (s). □ 



Assume that ip± is defined on the interval (—sq, so) and let $ t be the 
flow of Then define F±: (-s , s ) x R^ R 4 by 



(3.12) 



Ms) = ip±{s)s' 



(3.13) 



F ± (s,t) :=$ t ( 7± (s)). 



Writing this as 



F±{s, t) = w\(s, t)vx + w 2 (s, t)v 2 + w 3 (s, t)v 3 + w 4 (s, £)v 4 



FIRST ORDER NON-LINEAR PDE IN TWO VARIABLES 19 

from the discussion of differential equations for the characteristics we 
have 

( Wl (s,t) =se a9lt , w 2 (s,t) = ±se~ bd2t , 

w 3 (s,t) = (f±(s)e- ad3t , w A {s,t) = i)±{s)e mt ) 

where \0i — 1] < e < 1/2. We are now going use u = wi(s,t) and 
v = W2(s, t) as new variables and express ws(s, t) and u>4(s, t) in terms 
of these variables. As W2(s,t) = ±se~ W2t this calculation breaks up 
corresponding to the choice of signs. To simplify notation we will do 
the calculation in the first quadrant of the wv-plane, which corresponds 
to choosing + in the definition of W2, using the functions ip + , 
and restricting to s > 0. The calculations in the other quadrants are 
identical. In the first quadrant 

u = w^s, t) = se a9lt , v = w 2 (s, t) = se~ W2 \ 

with s, u, and v positive. Solving for s and t we get 

S = u b/(a+b) v a/(a +b ) e a b ( e2 - ei )t f = 1 ^ 

a6 1 + b6 2 V'. 

Thus 



t > u 



l/(ae 1 +W 2 ) /u\(l/(o+&))+/si 

where 



1 1_ (1 - eja + (1 - 6 2 )b 

Pl ~a9 1 + b9 2 a + b~ (a + b)(a9 l + b6 2 ) ' 
As we are assuming that |1— 0j| < £ < 1/2 this leads to 

\ Pl \ < =: C x e. 

a + b 



This implies 



where 



— U a+b V a+b I [UV ) a+b J 
b a 1 _ 1 \ 02 

= U a + b V a + b [uv ) , 



P2 ■= [ -^-T+Pi ) ab(6 2 -e 1 ). 
a + b 



As |0! -6 2 \< 2e, |pi| < 2e/(a + 6) < l/(a + &), 

P2 < — —r =■■ C 2 e. 
a + b 
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Now 

w 3 (s,t) = V± (s)e- a9st 

= <p±(u^v^ (uv- 1 )" 2 ) (u«- 1 )te* + * 1 ) ( - aft0 
= tp±(u^v^ {uv~ l ) p2 ) {uv~ l )^ b+p '^ 

where 

P3 := -7 Pia#3- 

a + o 

Using the estimates \p\\ < 2e/(a + b) and |1 — 9 3 \ < e < 1/2 (so that 
1/2 < # 3 < 3/2) this gives 

, ea 2e 3 4ae _ 

Ps < — —t + — ^"o a = — = : c s £ - 
a + b a + o 2 a + b 



Using the form of tp + given by equation ( |3.2|) we have 

w 3 = Tp + (u^v^ (uv~ i y 2 ) (u^v^b (uv~ 1 ) P2 ^ (uv~ 1 )^" Tb ' +P3 ^ 
(3.15) 

/ b a / | \ /jo \ hi — a a(l+l) 1 \ 7 _|_ 

= Tp + (u°-+ b v a + b [uv l y )u a + b v a + b (uv~ ) p2+p3 . 
Likewise 

w 4 (s,t)=^ + (s)e b9 * t 

= ^ + {u^v^ {uv- l ) p2 ) (uv- 1 )^^ 64 
= ip+^u^v^ (uv~ 1 ) P2 ) (uv~ 1 )("+ b ' +p ^ , 

where 

(04 - 1)6 

p 4 = -r- + PibOi. 

a + b 

Thus 

eb 2e 3, Abe 

P4 — 7 "I r— = =: G 4 £. 

a+o a + o 2 a+b 
Using the form of given by equation ( |3.12| ) we have 

wa = ip + (u'^v'^ (uv~ 1 ) P2 ) (u^v^ (uv~ 1 ) P2 ^ (uv~ 1 )^~ +p ^ 
(3.16) 

y b a , i v q<~, b(l + l) al — b ~\ \ J _l_ 

= %l) + (u a + b v a + b \UV ) )U a + b v a + b [uv~ ) tp2+p4 . 
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3.6. Definition. Let kbe a positive integer, m a positive real number, 
and U a connected neighborhood of the origin in R 2 . Then we denote 
the class of functions defined on U that are C k off of the coordinate 
axes and which vanish to order m along the axes by Caxes ^ — Caxcs ^ ( C/) • 
Explicitly h G Cix^\u) iff h is C fc on U s {m> = 0} and for any 
compact subset V of {/ there is a constant C = Cy so that 

\h(u,v)\ < C\u\ m and \h(u, v)\ < C\v\ m 

hold for all (u, v) G V. □ 

Using this definition we now summarize and extend the calculations 
above. 



3.7. Proposition. Let H: R 4 — > R o/ c/ass C°° 6e gwen fry pTi 
Ze£ </?-t : R — > R be C k (with 1 < k < oo) functions that vanish 
to order I at the origin (as in (\3.3[) ), and let ip±: R — > R and the 
curves j± : R — > R 4 6e jwen 6y Lemma \3. 5] . Let wi(s, t), . . . , ^4(5, t) 
be the solutions of the characteristic equations given by (\3.14 ) and de- 
fine u(s, t) — wi(s, t), v(s, t) = w 2 (s, t). As above, we can then solve for 
s and t in terms of u and v. Then define functions on a neighborhood 
of the origin in R 2 by 

g(u,v) := w 3 (s,t), h(u,v) = w 4 (s,t), foru,v^0, 

and extend this to the coordinate axes by 

g{u,0)=g{0,v) = h{u, 0) = h{0,v) = 0. 

Let m be any real number with 

{bl — a al — b 
-, - 
a + b a + 

Then there is a connected neighborhood U of the origin such that g,h G 

in a small neighborhood of s = 0, then g and 
h vanish in a neighborhood of the coordinate axes, and are C k . The 
surface 

N 2 = N 2 (ip + , (pJ) := {mvi + v\ 2 + g(u, v )v 3 + h(u, v )v 4 : (u, v) G U} 

is an invariant Lagrangian submanifold of R 4 for some sufficiently 
small neighborhood U of the origin in R 2 . 

Proof. Most of this follows at once from the calculation above. The 
flow $ t of the characteristic vector field is C°° and the curves 7± are 
C k . Therefore Wi(s,t) is a C k function of (s,t) for s ^ 0. So by the 
inverse function theorem s and t are C k functions of (u, v), and g and 
h are C k functions off of the coordinates axes u = and v = 0. Along 
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the u = axis the formula ( 3.15 ) for w 3 in terms of u and v implies 
that 

\g(u,v)\ = \w 3 \ < C\u\^ (uv' 1 ) ]lp2+p31 
on any compact set. We have 

\lp2 + Ps\< (lC 2 + C 3 )e. 

Since > m we can, by choosing 5 small enough, make e so small 
that — (IC 2 + C 3 )e > m. But then (for a possibly larger value of 
C) 

\g(u,v)\ = \w 3 \<C\u\ m . 
Similar considerations, using ( p,15|) and (|3.16|) , show 



\h(u,v)\<C\u\ m , \g{u,v)\ < C\v\ m , \h{u, v)\ < C\v\ m . 

While this only shows that this estimate holds in a small neighbor- 
hood of (0, 0), we can extend it along the u = follows. From 
the construction the subset iV 2 := {wvi + t>v 2 + g(u, v )v 3 + h(u, v )v 4 : 
(u, v) G V} is invariant under the local flow of $ t . Likewise the uv- 
plane, realized as R 2 ^ := {wvi + v v 2 : {u, v) G R 2 }, is invariant 
under the flow of $t, as is the line defined by « = 0, realized as 
R„ := {t>V2 : v G R}. The estimates above imply that in some small 
neighborhood of the origin iV 2 and R 2 v have contact to order m along 
Ti v . But the flow $ t is smooth and therefore preserves the order of con- 
tact. Thus invariance of iV 2 , R 2 v and Ti v under $t shows that A^ 2 and 
R 2 „ have order m contact along H v at all points that can be moved 
close to the origin by $ t . Similarly, we can extend the estimates along 
the v = axis. 

Finally, from the construction we have that iV 2 C {H = 0} and 
that the characteristic vector field is tangent to iV 2 . Let X be a 
vector tangent to iV 2 , then dH(X) = 0, because iV 2 C {H = 0} and 
therefore H is constant on iV 2 . From the definition of we have 
u(Zh,X) = -dH(X) = 0. As £ H is tangent to A^ 2 and T(N 2 ) is two 
dimensional, this implies that u = at points of A^ 2 where ^/0. 
But £h has an isolated zero at the origin and therefore by continuity 
^1^2 = m some neighborhood of the origin. This shows that A^ 2 
is an invariant Lagrangian surface near the origin and completes the 
proof. □ 

Recall that the coordinates x,y,p,q we have been working with are 
not the standard coordinates on R 4 , but rather the coordinates x, y, p, q 
that reduced H to the Sternberg normal form. We now need to trans- 



late Proposition 3.7 back into a result about solutions to the original 
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equation H(x,y, z x , z y ) = 0. However as the notion of being an in- 
variant Lagrangian surface is invariant under a symplectic change of 
coordinates, this translation does not involve much calculation. 

Let N 2 be the invariant Lagrangian surface from the conclusion of 
Proposition IIO. Then from the construction of N 2 the two curves 



ai and a 2 in R 4 defined by ( p.lGQ lie in N 2 . Thus a'^O) = vi and 
a' 2 (0) = V2 is a basis for T(N 2 )q. Then (recalling that vi = e% and 
v 2 = e4 at the origin) the condition ( |3.5| ) and the implicit function 
theorem imply that by restricting the size of N 2 we can assume that 
it projects over a neighborhood of (0,0). Again letting 7r: R 4 — > R 2 
be the natural projection, we then have that the two curves 7r o a x 
and 7r o a2 cross transversely at the origin. Let a\(t) = 7r(ai(i)) and 
o>2{t) = 7r(a2(i)). Then these curves only depend on the curves ai and 
a2 and thus only depend on the Sternberg normal coordinates x, y, p, q 
and the function H, but are independent of the invariant Lagrangian 
surface N 2 . As Si and a.2 are C°° and cross transversely at the origin, 
there are C°° local coordinates u, v centered at (0, 0) so that for t near 
zero 

w(a 2 (t))=0, £?(ai(t)) = 0. 

3.8. Definition. Let A; be a positive integer, m a positive real number, 
and U a connected neighborhood of the origin in R 2 contained in the 
domain of the coordinates u, v. Then denote by C^_^ xes (U) the class of 
functions defined exactly as in Definition |3l| but with the coordinates 
u and v replacing the coordinates u and v. □ 



With this definition we can now translate Proposition !T7 into a 
statement about solutions to the differential equation H(x,y, z x , z y ) = 
0. 



3.9. Theorem. Use the notation of Proposition ^/^ and assume that (\3.j\ 
and (\3.5) hold. Then given the C k (with 1 < k < oo) functions 



(p±: R — > R there is a solution z to H(x,y, z x , z y ) = defined in a 
neighborhood U of the origin that has N 2 (ip + , <£>_) as its jet near the 
origin. Moreover, for any real number m with 

bl — a al — 



m < mm 



a + b ' a + b 



we have z e c£ e + ^™ +1) ([/). 

Proof. That there is a solution z with iV 2 = N 2 (ip + , ip_) as jet follows 
from Proposition p72| . As iV 2 = {(x, y, z x , z y ) : (x,y) G U} there is 
a gain in regularity of one derivative in going from the jet iV 2 to the 
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function z. As N 2 has regularity Cixes , we see that z has regularity 

J u-v axes 



C-a-v aiTs • This completes the proof. □ 



3.10. Remark. Choosing <p + and to vanish to infinite order at the 
origin we see that under assumptions ( |3.4j ) and ( |3.5| ) any saddle type 
solution to H(x,y, z x , z y ) = will have an infinite dimensional family 
of smooth deformations. □ 



4. Formal power series and non-existence of smooth 

solutions 

In this section we will use power series methods to investigate the 
local existance and regularity of solutions to the eikonal equation, 

(4.1) 4 + z 2 y = h, 

near a point where h = 0. We assume that h is smooth and that it has 
a non-degenerate Hessian at the zero we are considering. We make, 
if necessary, an affine change of coordinates of the type a translation 
followed by a rotation, to bring the series expanion at this zero into 
the form 

(4.2) h(x,y) = a 2 x 2 + b 2 y 2 + £ h m , n x m y n . 

m+n>3 

for some a, b > 0. 

4.1. Proposition. If z is a C 3 solution to 
and if a ^ b, then 

(4.3) z(x, y) = l -(±ax 2 ± by 2 ) + 0((\x\ + \y\) 3 ). 

If a = b then still holds after making a change of coordinates 

x 

y _ ' 

for some angle of rotation £. 

Proof. Suppose that z has a series expansion z(x, y) = ax + j3y + 
\(Ax 2 + 2Bxy + Cy 2 ) + 0((\x\ + \y\Y). Then z 2 x + z 2 = a 2 + (3 2 + 0{\x\ + 
\y\). But h vanishes at (0,0) which, along with ( f4.1|) , implies that 
a — f3 — 0. Evaluating the second order terms of z 2 + z 2 and equating 
them with the second order terms of (^4.2| ), we get the equations 

(4.5) A 2 + B 2 = a 2 , B(A + C) = 0, B 2 + C 2 = b 2 . 

Suppose first that a ^ b, then we must have B = since otherwise 
the middle equation would imply that A = —C which, together with 



JJ) such thatz(0,0) = 0, 



(4.4) 



X 


i — > 


y 





cos £ sin £ 
- sin £ cos £ 
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the first and the last equations, contradict the assumption that a ^ b. 
It then follows immediately that A = ±a and C = ±6, and so ( |4.3| ) 
holds. This holds also if a = b and B = 0, i.e. A = ±a and C = ±a 
with any combination of signs. 

If a = b and B^O, then any solution of (|4.5|) is of the form 

A e = ±aVT^¥, B e = off, C e = TaVT^, 

for some 6 E [—1, 1],0 ^ 0. For a given 9, the function z = ^(A 9 x 2 + 
2B e xy + Cgy 2 ) + • • • is brought to the form (|4.3| ) by the change of 
variables ( f4.4|) , if £ is chosen such that tan£ = 6/(1 + y/1 — 9 2 ). □ 

As the convex and concave solutions are well understood (these are 
the ones corresponding to leading terms ±(ax 2 + by 2 )/2), we will focus 
on the saddle solutions. Thus assume that z has Taylor expansion 

(4.6) z(x,y) = ^(ax 2 - by 2 ) + z m^ m y n . 

m+n>3 

(Note that, if a = b, we work with the rotated coordinates.) Then 
z x = ax + ^2 mz m, n x Tn ' 1 y n 

m+n>3 

so that 

2 2 2, \ ^ ran,/ \ ^ m-1 

z x = a x + 2maz m ^ n x y + I mz m, n x y 

m+n>3 \m+n>3 

n >3 mz m,nX m l y n ) we get terms 
mz m)n x m - 1 y n kz K ix k - l y l = mkz m ^z k ±x m+k - 2 y nM = mkz m , n z k)l x p y q , 

where p + q = m + n + k + l — 2 > k + I + 1 as m + n > 3 and 
likewise p + g>m + n + las&; + />3. Thus there are polynomials 
Pm,n m (zij) in the variables z^j with i+j < m + n such that 




Putting this all together (and setting P m+n = when m + n = 3) we 
get 

z\ = a 2 x 2 + J2 ( 2m ^ z m,n + Piy <m+n \ Zij ))x m y n . 

m+n>3 

Likewise, 

Zy bx -\- ^ nz m ^ n x y , 

m+n>4 
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so a similar calculation gives 



zl 



& V+ E {-2nbz m , n + Q^l <m+n \z h] ))x m y n . 



m+n>3 



Plugging these equations into (|4.1| ) we get the recursion 



(4.7) 2(ma - nb)z m , n = h m>n - {P^ <m+n] (z hJ ) + Q [ ^ <m+n] (z hJ )) . 

4.2. Theorem. Let h be a smooth function so that the Taylor expan- 
sion ofh(x, y) near (0, 0) is of the form ( \4-%\) , where a and b are positive 
and linearly independent over the rational numbers. Then at the level 
of formal power series there is a unique saddle point solution z to 
with leading terms of the form ( U-(\ ) 

Proof. If a and b are linearly independent over the rationals, then ma — 
mb for all m and n. Thus the equation ( |4.7|) uniquely determines 
all the coefficients of z as required. □ 

4.3. Corollary. With the same hypothesis as in the last theorem, there 
is a smooth function z so that z 2 + z 2 — h vanishes to infinite order at 
the origin (0, 0). 

Proof. This follows from the well known theorem of Borel (cf. || and 
J7|, Thm. 1.2.6. p. 16]) that given any formal power series there is a 
smooth (i.e C°°) function having the given formal power series as its 
Taylor series. Thus if z is the formal power series solution given by 
the theorem, then let z be a smooth function that has z as its Taylor 
expansion. □ 

4.4. Theorem. Let a and b be positive real numbers that are linearly 
dependent over the rationals. Then there is a smooth function h of the 
form (U-Sj ) such that there is no solution z of the saddle point form 



to the equation (\4-A) - ^ n particular, ifm andn are positive integers with 
ma — nb = and m + n> 4 and c ^ 0, then for h = a 2 x 2 + b 2 y 2 + cx m y n 
the equation t \4-A ) w ^ n °t have any solution in formal power series (or 
any C k solution for k > m + n.) 

Proof. If a and b are positive and dependent over the rationals, then 
there are positive integers m and n with m + n > 4 so that ma — nb = 
0. Therefore, if there is a solution z to ( |4.1|) of the form (|4.6|) , the 
equation (|4.7|) implies 

h _ ( p[i+j<m+n]( \ _|_ nXi+i<m+n]( \\ 

u>m,n — m ,n K^'hJ I ' Vm,n K^hjJJ " 

But then, by adding cx m y m to h, c ^ 0, we get a new h of the form ( f4.2| ) 
so that there is no solution to ( |4.1| ) in formal power series. For a 
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specific example, let h = a 2 x 2 + b 2 y 2 and z = \{ax 2 — by 2 ). Then 
(z ) 2 + (z ) 2 = h and the equation z 2 + z 2 = h G + cx m y n will have no 
solution in formal power series, if ma — nb = and 0. □ 

A. Appendix: The stable submanifold theorem and the 
Sternberg normal form 

A.l. The stable submanifold theorem. The following is a state- 
ment of the Stable Submanifold Theorem adapted to our applications. 
If k > 1 and X is a C k vector field on a smooth n-dimensional man- 
ifold M, we denote by $ t the flow of X. That is, $> t is the locally 
defined one parameter group of diffeomorphisms of M that satisfy 
|$t(P) = X(P). As X is C k , the function (P,t) i-> $ t (P) is also 
C fc (c/. [0, p. 230] or H Thm 5 p. 86]). A point P G M is a Ziyper- 
foo/ic critical point of X iff X(P ) = and the linearization of X at 
Po, L : T(M)p — > T(M)p , has no eigenvalue with zero real part. This 
implies that Po is an isolated critical point of X. Letting, as usual, 
e tL = X]fclo(^)/^!' ^ ne stable subspace of X at Po is the linear 
subspace of T(M) Po defined by 

T+(M) Po := {v G T(M) Po : lim e tL v = 0}. 

t— >oo 

The unstable subspace is likewise defined by 

T_(M) Po := G T(M)p : lim e*S = 0}. 

t— >— oo 

If no eigenvalue of L has zero real part, then basic linear algebra implies 
the direct sum decomposition 

(A.l) T(M) Po = T_(M) Po © T + (M) Po . 

When all the eigenvalues of L are real and distinct (which is the case for 
most of the applications in this paper) then T + (M)p is the span of the 
eigenvectors corresponding to the negative eigenvalues and T_ (M) p is 
the span of the eigenvectors corresponding to the positive eigenvalues. 

A.l. Theorem (Stable Submanifold Theorem). If k > 1 and P is 

a hyperbolic critical point of the C k vector field X, then there is a 
connected open neighborhood U of Pq such that 

N + :={PeU : lim $ t (P) = P } 

t— >+oo 

(the stable submanifold of X at P ) and 

N„:={P EU : lim $ t (P) = P } 

oo 

( the unstable submanifold of X at Pq) are both embedded C k sub- 
manifolds ofU . The tangent space to N± at P is T±(M) Po so dim N = 
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dim T±(M)p and therefore ( A.l\) implies that N + and AL intersect 



transversely and that dimiV + + dimAL = dimM. □ 

This was originally proven in various degrees of generality by Hadamard 
H , Liapunov and Perron |12|] . Modern presentations of the proof 
can be found in [H and flTJL Chap. 6]. 



A. 2. The Sternberg normal form. Let if be a smooth function 
on R 4 and £h the characteristic vector field of H as given by ( [2.1|) . 
Let P be a critical point of H (and thus a rest point of £h) and let 
L be the linearization of at P. Then the matrix of L is given by 
equation (|2.2j ). If the eigenvalues of L are real and nonzero, then the 
formulas ( |2.3| ) imply that they are of the from a, b, —a, —b, for some 
a, b > 0. With this in mind, we can give Sternberg's normal form for a 
Hamiltonian near a critical point. 



A. 2. Theorem (Sternberg |18[). Let H be a smooth function on R 4 
and P a critical point of H with H(P) = 0. Assume that the eigen- 
values of the linearization L of at P are a, b, —a, —b where a, b > 
and that 

(A. 2) a and b are linearly independent over the rational numbers. 

Then there are local coordinates x, y, p, q centered at P and a smooth 
function f{u,v) of two variables so that 

1. x,y,p,q are symplectic coordinates. That is, if uj is the symplectic 
form then 

uo = dp A dx + dq A dy . 

2. The function f satisfies f(u,v) = u + v + 0(u 2 + v 2 ) or more 
formally: 

/(0,0) = 0, /„(0,0) =/«((), 0) = 1, 

3. In these coordinates H is given by 

H= 1 -f(f-a 2 x 2 ,f-b 2 f). 



This follows from Thm. 9 in Sternberg [FSL p. 603] applied to the 



transformations given by the flow of a Hamiltonian vector field, using 
the considerations from Section 7 in Sternberg [|l7], p. 818 - 819]. 
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